Relativistic quantum walks 
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By pursuing the deep relation between the one-dimensional Dirac equation and quantum walks, 
the physical role of quantum interference in the latter is explained. It is shown that the time 
evolution of the probability density of a quantum walker, initially localized on a lattice, is directly 
analogous to relativistic wave-packet spreading. Analytic wave-packet solutions reveal a striking 
connection between the discrete and continuous-time quantum walks. 
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The "quantum random walk," first coined by 
Aharonov et al. [lj, is a quantum generalization of the 
classical random walk. Consider a walker moving on a 
one-dimensional lattice, taking steps left or right based 
on the state of a coin. Classically, if the coin is flipped 
after each step, this generates a diffusive random walk. 
If the coin is quantum mechanical, however, it can be 
put into a superposition, and rotated by applying a fixed 
unitary operator. Aharonov et al. showed that this quan- 
tum procedure (or algorithm) can generate displacements 
that, on average, are much greater than the classical ran- 
dom walk. 

This discrete-time quantum walk (DTQW) has been 
re-discovered and extensively analyzed in the context of 
quantum computation . Two key properties are the fol- 
lowing: (i) the standard deviation of the walker's position 
grows linearly in time ((Ax)t ~ t), in clear distinction to 
the classical random walk ((Ax)t ~ t 1 / 2 ), (ii) for proper 
initial conditions the walker spreads out symmetrically, 
with a nearly constant probability distribution save for 
two peaks located at x+ = ±ct (where c = l/\/2 for 
the "Hadamard walk" Q ) , beyond which the probability 
quickly goes to zero. 

An entirely different approach to "quantizing" random 
walks was initiated by Farhi and Gutmann |3(. Begin- 
ning with the differential equation for diffusion on a lat- 
tice, they performed an analytic continuation to yield a 
Schrodinger equation with a finite-difference Laplacian 
operator. This continuous-time quantum walk (CTQW) 
was used by Childs et al. to construct a special 
search algorithm that is exponentially faster than classi- 
cal methods. Other local search algorithms (with square- 
root speedup) have been studied using both the discrete 
and continuous-time quantum walks, often with similar 
results |(J. However, to this author's knowledge, no 
physical explanation has been proposed to explain the 
similar performance of these two quantum walks. 

Before connecting these two walks, recall the connec- 
tion between the DTQW and the Dirac equation. As 
discussed by Meyer Q, this goes back to Feynman's 



"checkerboard," a discrete space-time path integral that, 
in the continuum limit, generates the propagator for the 
Dirac equation in one dimension This is best seen 
in the following unitary representation |9(, in which the 
DTQW is written as the discrete mapping: 



tp R (n,T+ 1) 
4> L (n,T + 1) 



= U 



(1) 



tpR and ipL are wave functions on an infinite lattice, and 
U is the product of a conditional translation operator 
and a spin rotation 
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U= [~(I + a z )D + -(I-a z )D- 



(2) 



Here the Pauli matrices {/, o~ x , a z } act on the spinor com- 
ponents, and the translation operator D acts on wave 
functions as (Dip)(n) = ip(n — 1). The continuum limit 
is found by letting the position x = ne, D = e~ lep (p is 
the momentum), 6 — me (m is the mass), and the time 
t = er. Using the Trotter formula, the limit e — > (with 
p, m, and t finite) yields 



it/e 



iHnt 



(3) 
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where Hp is the Hamiltonian for the one-dimensional 
Dirac equation (with K= c = 1, p = —id x ) [Toj |: 

idt^(x,t) = H D ^(x,t) = (-ia z d x +a x m)y(x,t). (4) 

While quite elegant, the properties of this continuum 
limit have been largely ignored in the extensive analysis 
of the Hadamard walk (in which e~ l7 " Ty ^a z is used 
in place of e~ %e<Tx in J5J). The closest related work is the 
continuum limit of the Hadamard walk recently found 
by Knight et al. H2, but this and the corresponding 
Airy function solutions are significantly different from the 
Dirac equation. Another notable work is that of Meyer 
who studied some of the wavelike properties of quan- 
tum cellular automata, but not the uniquely relativistic 
properties explored here. Understanding these proper- 
ties may have importance for quantum algorithms; it has 
already been shown that massless Dirac operators can 
improve a continuous-time search algorithm 
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Here I use explicit solutions of (@J to illustrate that 
the quantum-walk probability distribution is analogous 
to the spreading of a relativistic particle. The term 
"relativistic" is taken to mean any evolution of a par- 
ticle with a maximum speed limit. The same character- 
istic spreading — both relativistic and nonrelativistic — is 
found from a new solution to the quantum walk equa- 
tions without going to a continuum limit. Finally, 
this solution is found to be analytically related to the 
continuous-time quantum walk, providing a new link be- 
tween these two relativistic quantum walks. 

First, it is important to note that, using the Heisen- 
berg equations of motion, wave-packet spreading for any 
dispersion relation u>(p) can be written as 



(Ax) 2 = (Axf + {Av) 2 t 2 , 



(5) 



where (A«)o is the standard deviation of the group ve- 
locity v(p) = duj(p)/dp [l4(. For the Dirac equation 
the dispersion relation is lo(p) = \J p 2 + to 2 , and thus 
v(p) = p(p 2 + to 2 ) -1 / 2 < c = 1, i.e. there is a maximum 
group velocity, which is of course the speed of light. 

While this linear quantum spreading (Ax)t ~ t is uni- 
versal, the presence of peaks of the probability distribu- 
tion (at x± = ±ct) depends on the initial localization. 
To show this, I construct an explicit time-dependent so- 
lution of Q by the following Fourier representation: 



j\[ r°° i / i 



ipx — (a+it)uj(p) 



(6) 

The prefactor P+{p) = I + Hu/uj(p) projects the spinor 
onto the positive-energy eigenstates of Hp , while the pa- 
rameter a in the exponential allows arbitrary localization 
in position. The integrals can be done analytically to 
yield 

m( t) - mAf ( s ~ lK i{ m s)[a + i(t + x)] + K (ms) \ 
[X ' ' ~ ^ V«- x ifi(m«)[a + *(t-a;)] +Jfo(ms) ) ' 

(7) 

where s — [x 2 + (a + it) 2 ] 1 ' 2 , the normalization factor is 



M = y/n/2m [Ki (2ma) + K Q (2rao)] 



-1/2 



(8) 



and K n is the modified Bessel function of order n • 

The probability density is shown in Fig. ^ for two val- 
ues of a at t — and t = 50. The nonrelativistic wave 
packet (with large a) spreads as a Gaussian, while the 
relativistic wave packet (with small a) spreads near the 
light-cone (x± — ±ct) at the speed of light. 

Another exact solution of the one-dimensional Dirac 
equation was found many years ago . These examples 
demonstrate the existence of positive-energy states of a 
relativistic particle localized beneath its Compton wave- 
length fll) j despite well-known claims to the contrary. 
Such states appear to require entanglement between the 
spatial and spinor degrees of freedom |l8j — this is shown 
for (7J below. 
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FIG. 1: (a) Non-relativistic (a = 5) and (b) relativistic (a = 
0.5) solutions of the one-dimensional Dirac equation. The 
probability density p(x,t) = (x, t)^(x, t) is shown at time 
t = (dotted lines, arbitrary units) and at t = 50 (solid lines). 
Other parameters are the mass m = 1 and the speed of light 
c = 1. 



While the resemblance between the probability distri- 
bution of the Hadamard walk (see, e.g. 0) and the rela- 
tivistic wave packet in Fig. 1 is quite strong, it begs the 
question: what about the non-relativistic case? For the 
quantum walk, this requires an initial superposition over 
the lattice. A larger spread in the position leads to slower 
spreading, as expected by the uncertainty principle. 

An analytic solution to the walk equations Q and J2J, 
covering both the relativistic and nonrelativistic limits, 
can be found using a similar procedure as above. I use 
the Fourier analysis of and J5J and let D = e~ lk , in 
which case 



U 



e lk cos ( 



— le sint 



suit 



cosf 



(9) 



This matrix has eigenvalues e ±lul ^ k \ where uj(k) satisfies 
the dispersion relation ^| 

cosu(k) = cos#cosfc. (10) 

The wave packet corresponding to ©-0 is 



NT 1/1 



ikn—{a.-\-ir)uj{k) 



( n ) 

The prefactor P+(k) = (e"^ fc ) — U) projects the spinor 
onto the "positive-energy" eigenstates of U, while the pa- 
rameter a in the exponential allows arbitrary localization 
on the lattice. This solution can be written as 

if \- — ( In{T-l-ia)-e- ie I n . x {T-ia) 
nn,r) - ^\ K j n{T _ 1 _ ia) _ e -ie In+1 ( T _ ia) 

(12) 

where the function I n (z) is defined by 

i r 

I n (z) = — I dk exp(ikn — iuj(k)z), (13) 
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DTQW: 6 = 3%/7, a = 22, x = 225 



0.03 




DTQW: 9 = 3ti/7, a = 2.2, x = 225 



S 0.02 

CO 

o 0.01 




0.8 



CD 

r o.6 

c 

CD 

E 

_CD 

°> 0.4 

CO 



0.2 



4 6 
a = a / tan(6) 



10 



FIG. 2: (a) Non-relativistic (a = 22) and (b) relativis- 
tic (q = 2.2) solutions of the one-dimensional discrete-time 
quantum walk (DTQW). The probability density p(n, r) = 
ip^ (n, r)ip(n, t) is shown at time r — (dotted lines, arbitrary 
units) and at r = 225 (dots), along with the Bessel-function 
approximation I15H (solid lines). Other parameters are the 
rotation angle 9 = 3tv/7, the mass m = tant9 ~ 4.38, and the 
speed of light c = cos 9 ~ 0.22. To compare with the Dirac 
solution, the parameters were chosen such that cr ~ 50 and 
a/m — a/ tan 8 ~ a. 



with w(fc) given by and the normalization factor is 



N = [2I o (-i2a)-e i0 h(-l-i2a)- 



'Ml- 



■i2a)]-V 2 . 
(14) 

At this point, a crucial approximation can be made: if 
cos 9 is small, replace oj(k) by its lowest order expansion 
from l|10f) : ~ 7r/2 — cos t9 cos A:. This replacement 

conveniently yields the same maximum group velocity 
(cos 9) and allows the following approximation to (|13|l : 



I n (z) ~ e -(»-^/ 2 J n (^ 



cos( 



(15) 



where J„ is the Bessel function of order n . Using this 
approximation in 1|12|) . as shown in Fig. [21 compares quite 
favorably to a numerical calculation of JIJ and J2J). The 
analogy between this and Fig.QJis remarkable, taking the 
"speed of light" for the quantum walk as cos 9. 

A few comments are in order. First, I have shown 
that evolution on the line for the DTQW, in the rel- 
ativistic case, has fronts that propagate at the maxi- 
mum speed c = cos 9, in close analogy to a solution 
of the Dirac equation. Heuristically, the criterion for 
a relativistic walker is for the the initial localization 
(Ax)o to be less than the effective Compton wavelength 
A = l/(mc) = l/sin(9, where the effective mass is given 
by m — [d 2 uj(k)/dk 2 ]^ = tan 0. To approximate the 
Hadamard walk, the appropriate choice is 9 — tt/4, lead- 
ing to c = l/v2> m — 1, and A = y/2 > 1. Thus, the 
initial condition most widely studied |Tll | , with the walker 
localized at one position, has (Ax)o ~ 1 < A, leading to 
a relativistic quantum walk. 



FIG. 3: The entanglement, in ebits, of the Dirac solution 
(solid), and the discrete-time quantum walk with 9 — 3tt/7 
(dashed), as a function of the scaling parameter a = a/ tant9. 
The entanglement measure is the spinor entropy using a 
base-2 logarithm. 



Second, the analogy between the wave packets of (J7J) 
and (|12|l extends beyond the probability distribution to 
the entanglement between the spinor and spatial degrees 
of freedom ^^|. The entanglement as a function of the 
initial localization is shown in Fig. [3] By including only 
positive-frequency terms in the wave function, the en- 
tanglement remains constant in time. Note that here, as 
in Fig. |2 I have used the correspondence between the 
localization parameters a — a/ taxi 0, found by compar- 
ing the dispersion relations near k = p = 0. As discussed 
above, highly localized positive-energy states become sig- 
nificantly entangled in the limit a — > 0. 

Finally, I note that the particular choice of the spin- 
rotation (coin) of the quantum walk analyzed above has 
simplified the calculation. As an example, the Hadamard 
walk's dispersion relation sinw (k) = sinfc/V2 G3 does 
not satisfy uj{k) ~ w(Q) + fc 2 /(2m) for small k, but rather 
uj(k) ~ fc/%/2 — fc 3 /(12V2). This expansion, the essential 
approximation used in [l2|. does not lead to an obvious 
nonrclativistic limit to the Hadamard walk. 

There is, however, both types of propagation — 
relativistic and nonrelativistic — for the CTQW, defined 
by 

id t ip(n, t) = -7 (ip(n - 1, t) - 2ip(n, t) + i{j{n + 1, t)) . 
An exact solution for this walk can be found as above 



i/)(n, t) = N{2tz)- 1 dfc e 4fc "-( Q + 4 *M fc ) 
= Ne- 2 i ia+it h n J n (2~/(t - ia)), 



(17) 



with the dispersion relation ui(k) = 27(1 — cos k) and nor- 
malization factor N = e 2la [Jo (—M'ya)] _1 / 2 . The rela- 
tivistic and nonrelativistic evolution for this case is shown 
in Fig. 4. This solution is strikingly similar to (|12|) [using 
the relation l|15|) ]. both visually and analytically, assum- 
ing equal maximum speeds c = 27 = cos 9. 
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CTQW: y = cos(3ji/7)/2, a = 22, t = 225 
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FIG. 4: (a) Non-relativistic (a = 22) and (b) relativistic 
(a = 2.2) solutions of the one-dimensional continuous-time 
quantum walk (CTQW). The probability density p(n,t) — 
ip^ (n,t)ip(n,t) is shown at time t — (dotted lines, arbitrary 
units) and at t = 225 (solid lines). The remaining parameter 
is c = 27 = cos(3vr/7) ~ 0.22. 

A Bessel-function approximation to the DTQW simi- 
lar to l|llfl -l|15l). was recently found by an entirely differ- 
ent method [Rj . The physical content of this approxi- 
mation, however, is revealed by the wave-packet analysis 
presented here: when cos 9 <Cl(ff~ tt/2), the dispersion 
relations of both the CTQW and DTQW have the com- 
mon form u)(k) — w(0) + c(l — cosfc), with the "relativis- 



tic" property of a maximum speed (v(k) = duj/dk < c). 
This equivalence is quite unexpected, since it is the 9 <C 1 
(cos6> ~ 1) limit of ((TDJ}, w{k) ~ Vk 2 + 9 2 , that leads to 
the Dirac equation. 

Despite this quantitative equivalence, there still ap- 
pears to be two qualitatively distinct approaches to quan- 
tizing a random walk. A possible resolution is simply 
to consider the CTQW as the discretization of the one- 
dimensional nonrelativistic Schrodinger equation, and 
the DTQW as the discretization of the one-dimensional 
Dirac equation. 

The Schrodinger equation can be considered the quan- 
tization (by analytic continuation) of the diffusion equa- 
tion for Brownian motion. The Dirac equation can also 
be considered the quantization (by analytic continuation) 
of the two-velocity model for the telegrapher's equation 
pof. This model describes a particle that moves with 
a constant velocity left or right, switching its velocity 
randomly at some constant rate (2l| . This latter pro- 
cess corresponds precisely to the coined classical random 
walk originally described above. From this point of view, 
the two quantum walks are not two different quantiza- 
tion methods, but rather equivalent quantizations of two 
different stochastic processes: one described by the diffu- 
sion equation and the other by the telegrapher's equation. 
That both lead to propagation with a maximum speed 
is the surprising yet simple consequence of discretizing 
equations on a lattice. 

I gratefully acknowledge P. R. Johnson and A. J. Dragt 
for key discussions at the beginning of this work. 
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